A New Object- Weighted Measure of the SmaU-Scale 
Velocity Dispersion 



Jonathan E. Baker and Marc Davis 

Astronomy Department, University of California, Berkeley, CA 94720 

G^ . 

Q>^ , Abstract. We describe a new statistic for measuring the small-scale 

velocity dispersion of galaxies directly from redshift surveys. This statis- 
pi ^, tic is based on the object- weighted statistic proposed by Davis, Miller, 

O ' & White (1997). Compared with the traditional pair- weighted velocity 

C/^ , dispersion, our statistic is less sensitive to the presence or absence of rare, 

' rich clusters of galaxies. This measure of the thermal energy of the galaxy 

distribution is ideally suited for use with a filtered version of the cosmic 
energy equation. We discuss the application of the statistic to the Las 
Campanas Redshift Survey. The low observed dispersion strongly favors 
' cosmological models with low matter density, Qm ~ 0.2. 
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1. Introduction 

Redshift surveys provide an accurate measure of the Hubble recession velocity 
Hor plus radial peculiar velocity (relative to the velocity vq of the observer): 



CZ = Hor + (Vpec - Vo) • f (1) 

' for large samples of galaxies. In the gravitational instability paradigm for struc- 

ture formation, peculiar velocities grow in response to the total (not just visible) 
amount of clustered mass. Measurements of the magnitude of these peculiar ve- 
locities are an important cosmological probe, with the potential to discriminate 
I cosmological models, to constrain the bias of the galaxy distribution, and to 

constrain models for structure formation. 

On relatively large (few Mpc and up) scales, peculiar velocities define a 
smooth flow field, the divergence of which is simply given by the galaxy over- 
density times the parameter /? Q^/b, where b = 6g/6m is the galaxy bias. 
A number of methods have been devised for combining redshift and peculiar 
velocity surveys to obtain /? from these large-scale flows (Strauss, this volume; 
Willick, this volume). These methods have not fully converged to a consistent 
solution; while POTENT prefers /3 ~ 1 (e.g., Sigad et al. 1998), other methods 
tend to give low values /3 ^ 0.5 (e.g., Willick & Strauss 1998). 

Turning to small scales of order 1 Mpc, we find that peculiar velocities 
are essentially thermal or incoherent. The kinetic energy contained in these 
motions can be used with a filtered version of the cosmic energy equation to 
yield an estimate of Qm/b'^i which is approximately equal to 

It has been known for some time that the thermal energy of the galaxy 
distribution is quite low, or equivalently the cosmic Mach number is rather high, 
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relative to our theoretical expectations. A simple way to see this is to take an 
A^-body simulation and plot what it would look like in redshift space. In a box 
of length 5000 km s"^, one sees very prominent, long "fingers of god", which 
tend to wash out the collapsed filamentary structures typical of real redshift 
surveys. The flow in the vicinity of the Local Group has also been measured to 
be remarkably cold, with a dispersion of only 60 km s~^, and no blue-shifted 
galaxies are seen outside the Local Group (Schlegel, Davis, &: Summers 1994). 
Govcrnato et al. (1997) showed that this is extremely hard to reproduce even 
in low-f7„j A^-body simulations. 

The pair velocity dispersion {(T12) is the traditional measure of small-scale 
velocities. Redshift surveys provide a set {rp,7r} of projected and radial sepa- 
rations in redshift space; ai2 can be estimated directly from the redshift-space 
correlation function ^z{fp, tt), which is a convolution of the real-space correlation 
function ^(r) with the pair velocity distribution function f(v) (the rms of which 
is (T12). An exponential f(v) has been found to fit well and is also expected from 
theoretical considerations. 

The CJ12 statistic was first applied to the CfA redshift survey by Davis & 
Peebles (1983), who measured 340 km s~^ at lh~^ Mpc scales. The fact that 
this dispersion was much lower than that of the = 1 simulations of Davis 
et al. (1985) was the original motivation for bias in the galaxy distribution. 
Gelb & Bertschinger (1994) later showed that no normalization of the standard 
Cold Dark Matter (SCDM) model could simultaneously match both the observed 
correlation amplitude and velocity dispersion. However, many authors pointed 
out that the pair statistic was not robust (e.g., Mo, Jing, &; Borner 1993; Zurek 
et al. 1994; Somerville, Primack, &; Nolthenius 1997; Guzzo et al. 1997; Jing, 
Mo, & Borner 1998). In particular, its pair- wise weighting makes (T12 very 
sensitive to the presence or absence of rare, rich clusters in the survey volume; 
the treatment of a few objects can greatly affect the result. The mean streaming 
motions of galaxies were also shown to have a considerable effect on the measured 
value of the dispersion. 

In this work, we present the application of a new, object- weighted statistic 
to the Las Campanas Redshift Survey (LCRS; Shectman et al. 1996). The 
LCRS is the largest existing redshift survey which is nearly fully sampled. It 
contains 26,418 galaxies with a median redshift of approximately 30,000 km s~^. 
The survey consists of six 1?5 x 80° slices (three at northern declinations and 
three in the south), and contains about 30 clusters. 

Even with such a large and well-sampled survey as the LCRS, the mea- 
surement of the pair dispersion au has been fraught with controversy. Fourier 
techniques (Landy, Szalay, & Broadhurst 1998) gave a low dispersion consistent 
with the old CfA value, but Jing & Borner (1998) showed that accounting for 
mean streaming motions increased the value to 570 it 80 km s~^. Jing et al. 
(1998) found that simulations with ^ 0.2 could reproduce the observed 

pair dispersion, but they required a somewhat mysterious anti-bias in the galaxy 
distribution (i.e., higher mass-to-light ratios in dense regions). 

Given the difficulty of reliably estimating the pair dispersion, it is clearly of 
interest to develop alternative statistics. Our statistic, called ai, is a modified 
version of the object- weighted statistic developed by Davis, Miller, Sz White 
(1997, hereafter DMW). We note that other alternatives to au have been sug- 
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gested. Nolthenius & White (1987) proposed the mean dispersion of groups 
of galaxies; our statistic has the advantage that it does not require assigning 
galaxies to groups and averaging over the internal motions. The pair disper- 
sion can also be measured as a function of local density (Kepner, Summers, &; 
Strauss 1997; Strauss, Ostriker, & Cen 1998), but unlike ai, this must be com- 
puted in volume-limited samples. This work discusses the application of ai to 
the LCRS; a more complete description of our results may be found elsewhere 
(Baker, Davis, & Lin 1999). 

2. The cJi Statistic 

We start from the definition of the single-particle object- weighted dispersion 
proposed by DMW. These authors applied this statistic to the UGC and IRAS 
1.2-Jy redshift surveys, which are much smaller than the LCRS survey. For the 
UGC catalog, they measured ai = 130 it 15 km s~^, much colder than their 
= 1 A^-body simulation, even when the simulation velocities were artificially 
cooled by a factor of two. 

2.1. Definition of cti 

Around each survey galaxy i in redshift space, we place a cylinder of radius 
Vp and half-length vi, with the axis of symmetry along the redshift (radial) 
coordinate. All the galaxies which fall within the cylinder are considered to be 
neighbors of galaxy i. We typically take Vp = lh~^ Mpc and vi = 2500 km 
s~^. We construct a histogram of the number of neighbor galaxies Pi{Av) in 
bins of velocity separation Av. We subtract the background distribution Bi{Av) 
expected for an uncorrelated galaxy distribution, and then average over the Ng 
galaxies to obtain the final distribution: 

D{Av) = ^ E [^^(^^) - Bi{Av)] , (2) 

where the weight for galaxy i is denoted by wi. 

In the original DMW formulation, the weight for each galaxy was simply 
given by its total number of neighbors in excess of background: 

w-^ = Y.[nAv)-B,{Av)]. (3) 

At; 

It is this factor which gives the statistic its essential object weighting, in contrast 
to the traditional pair-weighted dispersion. However, galaxies with fewer neigh- 
bors than the background had to be deleted from consideration, which biased 
the statistic towards higher density (and thus hotter) objects. 

We avoid the bias inherent in the original statistic by assigning Wi = 1 
for galaxies with less than one excess neighbor. In order to combine the high- 
and low-density galaxies sensibly, it is necessary to tabulate their distributions 
separately. We subtract off the tails of the distributions (within 500 km of 
vi) and then normalize each distribution so that the sum J2av D{Av) is propor- 
tional to the number of galaxies included. This allows the high- and low-density 
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Figure 1. Object- weighted velocity distributions for the LCRS. We 
show distributions for galaxies in low-density (dashed) and high-density 
(dotted) regions, and the combined distribution (solid). The widths are 
(Ti = 99, 207, and 126 km s^^, respectively. 

distributions to be combined so that they are weighted according to the number 
of objects assigned to them (see Figure ||). 

We measure the width ai of the resulting distribution by performing a 
fit to a model 

M{Av) = S^* f * E. (4) 

Here is the two-point correlation function ^(r) averaged in the cylindrical bins, 
the velocity distribution function, and E is the error distribution 
of the redshift measurements (for the LCRS, a Gaussian of rms 67 km s-^). We 
find that the exponential form for / fits much better than a Gaussian for both 
the data and A^-body simulations. We have defined ai so that it is a measure of 
the one-dimensional dispersion of the motion of individual galaxies, with bulk 
flows filtered out on scales larger than lh~^ Mpc. We note that the rms of the 
distribution / is a\\f2^ which DMW distinguished from (j\ by calling it ai. 

2.2. Cosmic Energy Equation 

The cosmic energy, or Layzer-Irvine, equation is a differential relation between 
the kinetic and potential energies of the mass fluctuations in an expanding 
universe. It has been shown that for self-similar cosmological clustering, the 
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equation reduces to an algebraic expression: (Upec) ~ 9^mHQJ2, where J2 = 
/ r^(r) dr measures the potential energy and g ~ 0.3. 

The difficulty in applying this equation to measure the mass density arises 
because the kinetic energy term and J2 are very poorly constrained on large 
scales. We therefore consider a filtered version of the equation, including con- 
tributions only from small scales. The simulations of DMW showed that ai is 
a good measure of the small-scale kinetic energy and can be used to estimate 
^m/b'^ (the factor of 6^ arises because we can only measure J2 for the galaxies 
rather than the underlying mass). 

2.3. iV-body Simulations 

We have completed several cluster-normalized A^-body simulations for compar- 
ison with the LCRS. The cosmological parameters of our models are given by 
Baker et al. (1999). They include a standard = 1 (SCDM) model, a 
tilted variant with n = 0.8 (TCDM), a flat model with a cosmological constant 
J)a = 0.7 (LCDM), and an open model with = 0.3 (OCDM). 

The models were evolved using a P3M code (Brieu, Summers, Sz Ostriker 
1995) and a special-purpose GRAPE-3AF board (Okumura et al. 1993), which 
is hardwired to compute the Plummer force law very quickly. The board is 
attached to a Sun SPARC workstation which computes the long-range con- 
tributions to the forces. We arc able to complete one cosmological run with 
64^ particles on a 128^ mesh in approximately one CPU-day. Our box size is 
L = 50h~^ Mpc to match the length of the redshift-space cylinders used in the 
o"! analysis. 

We apply the same statistical procedure for computing cji for the LCRS 
to the particles in the simulations. Of course, galaxies may in general have 
a different velocity dispersion from the underlying mass, and the problem of 
identifying "galaxies" in the simulations is therefore an important one. 

We first apply the standard friends-of-friends algorithm for defining halos, 
with a linking length of 0.2 mesh cells. With our relatively poor mass resolution, 
this procedure leads to a serious and well-known over-merging problem, yielding 
halo correlation functions which are much too low on small scales. To remedy 
this situation, we subdivide large (N > Ng) halos by drawing individual particles 
from them at random, with a probability oc N~'^. For a > 0, this yields a mass- 
to-light ratio which increases at small scales; this is required because a = leads 
to correlation functions which are too steep in many models. 

We choose the parameters a and Ng to yield "galaxies" which match the 
correlation function and number density of the LCRS (about 2500 per simula- 
tion volume). We find that a ~ 0.25 is typically the best choice, and Ng ~ 80, 
corresponding to a mass of W^^Q,.nh~^ Mq. Although by selecting individual 
particles as galaxies we are including the internal dispersions of galaxies, these 
velocities are small compared to the dispersions of the large clusters which are 
split. We have also measured ai for galaxies drawn from a large Virgo simu- 
lation (Benson et al. 1999) using semi-analytic techniques, and we find results 
consistent with our simulations. 
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Figure 2. Velocity distributions for the LCRS data (histogram) and 
best-fitting model (solid curve) with ai = 126 km s~^. Also shown are 
a model with no broadening (cji = 0; dotted) and the distribution for 
the mass in our SCDM simulation (ai = 310 km s~^; dashed). 

3. Results 

Based on the mean of the six LCRS slices, we measure ai = 126 it 10 km s~^, 
where the quoted error is the standard deviation of the mean for the slices. The 
statistic is quite robust, with similar D{Av) distributions for each of the slices. 
The fit is quite good, with xt — 117/96 = 1.22, where we have estimated the 
errors from the standard deviation of the six slices. The distribution is plotted 
in Figure ^. 

The exponential model also provides an excellent fit to the D{Av) distribu- 
tions in the simulations. We find that the mass in the SCDM and TCDM models 
is much too hot, with dispersions of over 300 km s~^. The LCDM and OCDM 
models yield lower dispersions ~ 200 km s~^; the OCDM model is slightly hot- 
ter than the LCDM. The halos are somewhat cooler, with velocity biases in the 
range 0.7-0.9. The LCDM halos are the best match to the LCRS dispersion, 
with o"! = 143 km s~^. 

Combining the LCRS and A^-body halo results, we can solve the cosmic 
energy equation for (the bias factor drops out because we have chosen halos 
which match the LCRS correlation function). We obtain similar results from 
each of the four cosmological models, with in the range 0.15-0.25. If we 
combine the LCRS and A^-body mass, we obtain ^l^/b, approximately equal to 
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Figure 3. Dependence of a\ on the cylinder length and radius. Filled 
points with error bars are for the LCRS, open points and dashed line 
are for galaxies in our LCDM simulation. 



/?. The result is 0.3-0.4 for the two high-density models, and 0.4-0.6 for the two 
low-density models. 

As we increase the radius of the cylinders used to measure cti, we find 
an interesting discrepancy between the models and the LCRS data (Figure ^). 
The LCRS a\ decreases modestly as the scale is increased, but all of the models 
show an increase with scale. It is at present unclear whether this discrepancy is 
a numerical artifact of the simulations or a real physical effect. 



4. Conclusions 

We have applied a new object-weighted, unbiased measure of the small-scale 
velocity dispersion to the Las Campanas Redshift Survey. We derive a single- 
particle dispersion g\ = 126 ± 10 km s~^. Our statistic has considerable advan- 
tages over the traditional pair dispersion ai2] namely, it is less sensitive to rare, 
rich clusters of galaxies. Our statistic should play an important role in analyses 
of future redshift surveys. 

When compared with the LCRS data, cluster-normalized 0^ = 1 A'^-body 
simulations are far too hot on small scales. We find strong evidence for a low 
density VL ~ 0.2, and we derive consistent values of the density parameter from 
a number of different models; these results are described more fully elsewhere 
(Baker et al. 1999). In addition to constraining the mass density, the cJi statistic 
applied to upcoming surveys should provide important constraints on the galaxy 
bias and evolution of structure. 
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